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Abstract 

We present a class of anisotropic brane configurations which shows BKL oscillations near their 
cosmological singularities. Near horizon limits of these solutions represent Kasner space embedded 
in AdS background. Dynamical probe branes in these geometries inherit anisotropies from the 
background. Amusingly, for a probe M5 brane, we find that there exists a parameter region where 
three of its world- volume directions expand while the rest contract. 
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Introduction 



AdS/CFT correspondence proposes an equivalence between the string theory in AdS space and 
a conformal field theory on its boundary. One of the remarkable features of this correspondence 
is the fact that it relates a strongly coupled theory with a weakly coupled one. Consequently it 
provides us with a way to tame the non-perturbative region of one by performing computations 
on its dual. Due to strong gravitational fluctuations, physics around cosmological singularities is 
dominated by non-perturbative effects and one hopes that the AdS/CFT correspondence would 
shed some light into it. Indeed, in recent years, we have witnessed several important investigations 
where attempts were made to find the signatures of these singularities in their gauge theory 
duals. Expectation is that the dual theory evolution might be able to provide a sensible quantum 
description of these singularities. Successes have been varied, please see the references p] - [5]. 

Inspired by this line of developments, in this letter, we search for D brane solutions in ten 
dimensional type IIB theory where the world volume metric expands anisotropically and show 
instabilities within their supergravity descriptions. We find that appropriately tuning the five 
form field strength, it is possible to construct a D3 brane with four dimensional Kasner like world 
volume. Along with a time-like singularity at r = 0, the metric shows an additional cosmological 
singularity at t = 0. Perturbation around t = generates an analogue of Belinskii-Lifshitz- 
Khalatnikov (BKL) oscillations. The near horizon geometry of this brane reduces to that of a 
Kasner-universe in AdS space plus a five sphere along with an appropriate five form field strength. 
We further probe the geometry by a dynamical D3 brane whose world- volume inherits anisotropic 
expansion/contraction along with a BKL like oscillation. Similar solutions can be constructed 
even within eleven dimensional supergravity. As an illustrative example, we discuss the case of 
M5 brane. The near horizon geometry is now a six dimensional Kasner space. A dynamical 
probe M5 brane in this space-time again acquires an anisotropic expansion in some directions 
and contraction along some. Amusingly, we find that it is possible to tune parameters in such a 
manner that three directions expand and the rest contract. Close to the cosmological singularities 
supergravity descriptions of all these solutions are expected to break down. We hope that the 
gauge theory description would shed some light on the physics near the singularities. 



D3 brane with anisotropic time-dependent world volume 

Besides the static D brancs of odd space dimensions, the IIB string theory admits time dependent 
branes. Consider, for example, the case of D3 brane. The equations of motion following from the 
relevant part of standard IIB supergravity action 
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has the following forms: 
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The numbers a, /3, 7 can be organized in an increasing order a < j3 < 7 and they vary in the 
range 
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These numbers can also be parametrized as 

—u 



a(u) 



1 + u + u- 



1 + u 



1 + it + u- 



■j 7(«) 



U + 7T 

1 + it + n 2 ' 



(6) 



where the Lifshitz-Khalatnikov parameter u > 1. Further, values u < 1 lead to the same range 
as 

(7) 



a(-) = a(u), /3(i)=/?(u), = 7 («). 
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The five form charge can be calculated by integrating *i ? 5 over the transverse space and it turns 
out to be time independent. 

In our convention, the extremal D3 brane is represented bya = /3 = 7 = and is not 
continuously connected to the above solution. Unlike extremal D3 brane, this solution breaks all 
the supersymmetries of IIB theory due to its explicit time dependence. The Kretschmann scalar 
for the metric is given by 
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In writing the above equation, we have used the condition @. It has a time-like singularity at 
r = at any finite time. It, further, has a cosmological singularity at t = 0. 

In the large r limit, equations in ([3]) reduce to a four dimensional Kasner solution plus a flat 
six-dimensional part. Within the Bianchi classification of homogeneous spaces, the Kasner metric 
corresponds to choosing all three of the structure constants to be zero. A generic perturbation 
near the singularity breaks these constraints generating Belinskii-Lifshitz-Khalatnikov (BKL) 
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oscillations [6|. To briefly illustrate the BKL oscillation, appropriately generalized to our context, 
we replace the world volume metric on the brane by type IX homogeneous space. 



To this end, let us consider the brane configuration of the form 

ds 2 =(! + -) 2 [ - dt 2 + (a(t) 2 lilj + bitfniimj + c(t) 2 n i n j )dx 1 dx 3 
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with the anti-symmetric five form field and the scalar 



_ 4r 3 Z 4 a(t)6(t)c(t)sin(x) 

= 0. (11) 

Here Zj,mj,nj are frame vectors. For IX metric, all the three structure constants are 1 and the 
simplest choice for the frame vectors is 

li = (since, — cosz since, 0), = (cosx, sinz since, 0), rij = (0, cosx, 1). (12) 

The coordinates run through values in the ranges < x < 7T, < y < 2tt, < z < 4710. The 
above configuration ([91- Hip is a solution provided they satisfy IIB equations of motion (|2|). This 
requirement leads to the following differential equations for a, b and c. 
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[(c 2 - a 2 ) 2 - 6 4 ] , 

[(a 2 - 6 2 ) 2 - c 4 ] , 
a n bn c+t 

— + -T- + — = , (13) 
abc 

where the subscript indicates derivative with respect to t. These are exactly the equations 
responsible for generating standard BKL oscillations. Consequently, the brane world-volume 
metric will oscillate with negative powers of t oscillating from one direction to another. In the 
next paragraph, for the sake of completeness, we give a brief analysis of this oscillation. 

To proceed, first we notice that if all the expressions on the right hand side of (|13p are 
small in some region, the system will have a Kasner-like regime with 

a~t a ,6~t^, c~t 7 , (14) 



* In all our discussion, we will closely follow [7j. [9] also has a lucid review of BKL oscillations for types VIII 
and IX spaces. 

^For type I spaces, in which Kasner metric belongs, the right hand sides of all the equations in (|13|) would have 
been zero. This is due to the fact that all the structure constants are zero for type I spaces. 
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where a,/3,7 satisfy constraint as in However, now since a is negative, close to t = 0, a 4 
term in the right hand sides of (j 13f) will start dominating. It is useful to write these equation in 
terms of new variables defined as 

a = e p , b = e q , c = e s , e p+q+s dr = dt. (15) 

In the vicinity of t = 0, (fT3|) reduces to 

Ptt = ~\ eAp i Qtt = s TT = ^e 4p , (16) 

where the subscripts r indicates derivative with respect to r. The solution of these equations 
should describe the evolution of world-volume metric from the initial state of Kasner metric. In 
terms of the new variables, this is equivalent to 

p T = a, q T = P, s T = 7. (17) 

Note that the first equation in (|16p can be interpreted as a particle moving in the presence of an 
exponential wall-like potential. Due to the reflection from this barrier, particle will move with 
p T = —a. However we see from (fl~6l) that p T + q T and p T + s T are constants. So we get 

q T = f3 + 2a, s T = 7 + 2a. (18) 

These lead to 

e P = e ~ aT , e q = e ^ +2a)T , e T = e^ +2a)r and t ~ e (1+2a)r . (19) 
In terms of the original variables, we can re- write the above as 

\a\ ff-2|a| 7-2|q 

a = i , 6 = ti-2|«l, c = t 1 - 2 ! Q l. (20) 

Therefore the action of the perturbation results oscillations between one Kasner regime to another 
with negative power shifting from a to b to c, inducing BKL oscillations on the world-volume 
of the D3 brane. We should however note that near the curvature singularity, string action 
receives higher derivative gravitational corrections. Consequently, the nature of the singularity 
and behaviour of its perturbation may get substantially modified. Solutions may also get modified 
when one introduces other matter fields into the theory. If they are represented by a perfect fluids 
on the world-volume, with pressure and energy density related as p = ujp, it can be argued that 
for cj < 1, BKL oscillations still persist. However, situation changes drastically for u) = 1, namely 
for the stiff-matter (a massless scalar field for example). A general discussion on these issues can 
be found in [U [TO]. Indeed it is easy to check that in our previous solution, one can introduce a 
dilaton with a profile 

= Alogi. (21) 

The metric and the form field remains same as before. However, the exponents now satisfy new 
constraint relations : a + /3 + 7 = 1, a 2 + /3 2 +7 2 = l — A 2 . The changes in these relations allow 
BKL oscillation for a finite time and the system finally reaches an attractor. 

We now proceed to study the metric in the near horizon limit r — > 0. In this limit, the 
metric reduces to 
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with 

4tr 3 tr 4 
Ftxyzr = giving potential C txyz = (23) 



The Kretschmann scalar is, 



We call it a Kasner-AdS space. This Kasner-AdS solution separately satisfies five dimensional 
Einstein equation in the presence of a negative cosmological constant and was found in [TTJ in 
the context of brane-cosmology. 



Probing with a D3 brane 

In this section, we will probe the geometry ([3]) with a D3 brane. Distance of the probe brane 
from the source now behaves as a scalar whose explicit time dependence can be determined via 
a dynamical equation. We take the world- volume directions of the D3 brane as £ = (t,x,y,z). 
The world-volume action of the D3 brane in background geometry, ([3]) takes the form : 

S = T J d*Zy/-det G a p+T J d*f <7 4 , (25) 

Here G a /s is the induced metric on the world-volume and C4 is the pull-back of the background 
4-form potential. T is the brane tension. We turn off all other fields on the brane. 
The Lagrangian can be cast in the form : 



L = y/A(t,r) - B(t,r)r 2 - C(r), (26) 

where, 
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The equation of motion for r(t) is the Euler-Lagrange equation derived from (I27P : 

r 4 (l 4 + r 4 f[tr + r] - (Z 4 + r 4 ) 3 r 3 - 2tZ 4 r 3 [3(Z 4 + r 4 )r 2 - 2r 4 ] 

-4t/ 4 r [r 4 - (I 4 + r 4 ) r 2 ] * = (28) 



Here dot represents derivative with respect to t. Once this equation is solved with appropriate 
boundary conditions, the metric on the probe brane is uniquely determined. We will carry out this 
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computation in this section. However, owing to the explicit time dependence in the background 
geometry, we find that the dynamical equation can not be solved analytically. Fortunately, it is 
not hard to find numerical solution and a typical behaviour is shown in the figure 1. 
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Figure 1: Plot of r as a function of time, t. 



The functions that govern the anisotropic expansions in three spatial directions are t a f(r), 
lP /(r), i 7 f(r), where 

/W=(l + £) 4 " (29) 

In order that the near horizon geometry is an AdS, as mentioned earlier, a,/3,7 must satisfy the 
constraint : a 2 + /3 2 + 7 2 = a + (3 + 7 = 1. This means, once we specify one of the three, say, a, 
the other two are automatically fixed : 

P = - (l - a + \/-3a 2 + 2a + l) 

7 = -(l-a - \/-3a 2 + 2a + l) (30) 



Ideally, in cosmology, one defines cosmological time, r\ with which the metric on the probe 
brane takes the form : 

dSLne = -drf + (l + ^) 2 (t 2a (ri)dx 2 + t 2 ^r,)dy 2 + t^(r,)dz 2 ). (31) 

with 

The behaviour of time, t as a function of r/ is depicted in figure 2. 
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Figure 2: Plot of t as a function of rj. 

At this, we plot the functions t a f(r), t* /(r), t" 1 f(r) as functions of r\ parametrically. 
Here f(r) is defined through f)29f) . One can tune the values of a, /3, 7 consistent with the 
Kasner constraints so that one of them goes down to zero (decelerating) while two of them go up 
(accelerating) with cosmic time and vice versa. One such plot is given in figure 3. 
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Figure 3: The functions, /1 = t a {r]) f(r), fa = t^{rj) /(r), /3 = t~<(r}) f(r), with f(r) given in 
(|29j) are plotted as functions of r\. a, (3 and 7 are 0.7, .632, -0.332 respectively The plot of f\ is 
in red, and that of f2 and are in green and black respectively. 
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The M5 Brane 



Our previous discussion can easily be extended to eleven dimensions. Here we discuss the case of 
a M5 brane. We start with d = 11 supergravity action 

^■^/^(""i")' (33) 

which is a generic action for the bosonic part of d = 11 supergravity so long as we concentrate 
on static, flat translationally invariant p-brane solutions. 

The equations of motion arising from (I33j) admits a solution of the form : 



ds 7 



1 + - 



-dt 2 + t 2a Ux 2 + t 2a *dx 2 + t 2a Ux 2 + t 2 <dx\ + t 2 <dx\ 



1 + [dr 2 + r 2 dn 2 ], 



along with 
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provided a\ + «2 + + 04 + CK5 = 1 and a\ + a\ + a\ + a\ + a\ = 1. 

In the near horizon limit, i.e. r — >• 0, the metric and the non-zero component of the form 
field reduce to the forms : 
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and hence the potential is given by Ct Xl x2x 3 x4x 5 

We now make the following change of coordinates : 

w 2 = I. (37) 
With this, the metric in (1361) takes the form : 



ds 2 = g + p a i d -2 + pa * d -2 2 + pal d -2 + fa^-2 + ? a§^ + 4 l 2 ^ + Z 2 ^, (38) 

where Xj and t are suitably scaled versions of the coordinates, Xj and t respectively. It is worth 
mentioning in this regard that the scaling of the coordinates will not be the same because of the 
presence of different powers of t in front of dxf. This is a consequence of anisotropy. 
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Following our nomenclature, (j38|) is a metric of seven dimensional Kasner-Ads space plus 
a four sphere. For on = for i = 1, ...,5, this reduces to our known AdS7(2L) x S 4 (L) solution. 



Probing with a M5 brane 

In the same spirit as we considered the case of probe D3 brane, we now consider a probe M5 
brane in the background (f34"j) and ([35]) . 

In PST formalism [12], the world- volume action of M5 brane is given in terms of a gauge 
invariant 3-form field strength, = dA^ +C^\ where A^ is world- volume 2-form and 
target space 3-form. The world volume action in this formalism is written as : 



Sm5 = Tm5 J d 6 £ [£-DBI + £-KE + £-W z] , (39) 



where 



C-dbi = y —det [Gij +T~Lijj is the Dirac-Born- Infeld Lagrangian, 

Cke = — n^ klmn ^imnHjkpG pq diad q a is the kinetic piece for the 3-form, 

24 (da) 

r ijklmn 

o! 



'-ijklmn ±u ^l mn 



is the Wess-Zumino term. 

(40) 



Here S a @ is the induced metric on the world-volume, and are the pull-backs of the 3-form 
and 6-form background potentials respectively. T~L is defined as 

"a" is an auxiliary scalar field introduced in PST formalism to maintain manifest covariance. 



If we now take the world-volume directions of the M5 brane as £ = (t, xi, X2, £3, #4, x$), 
it can be explicitly checked that, in this "static gauge", there will be no component of in 
world-volume directions. We further simplify the system by turning off the world-volume 2-form, 
A^ . With all these taken into account, the full Lagrangian takes the simple form : 

L = y/A(t,r) - B(t,r) f 2 - C(r), (42) 
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where, 



A(t,r) 


= t 2 (l + 


B(t,r) 


= t 2 (l + 


C{r) 


tl 3 


Z 3 + r c 



(43) 

Here dot represents derivative with respect to t. The Euler Lagrange equation for r(t) is : 

2r 4 (Z 3 + r 3 ) 2 [tr + r] - 2r (Z 3 + r 3 ) 3 r 3 - 3t/ 3 r 3 [3 (Z 3 + r 3 ) r 2 - 2r 3 ] 

-6tl 3 A [r 3 - (Z 3 + r 3 )r 2 ] ' = . (44) 



In order to draw a cosmological interpretation of the solutions we obtain from ([44]), as 
usual, we go to the "cosmic time" coordinate, ij, in which the metric on the brane assumes a form 



z 3 \-* ' 5 



ds 2 brane = -d V (t) 2 + ^1 + — ) ( ^2t^( V )dx 2 ) , (45) 

with 



J=i 



d 4 = Jh + K 



dt V V t 



The functions that govern the expansion of the universe in the spatial world- volume directions of 
the brane are in this case t ai f(r), where 



/(»■)= U + • ( 47 ) 



We can choose a^s so that three of them are the same and mimics isotropic expansion in 
three directions. The other two are anisotropic. Such a situation can be parametrized as : 

ol\ = a2 = a% = p 
a 4 = r (V-15p 2 + 6p+l-3p + l 

a 5 = 2 (~V-^P 2 + 6p + 1 - Sp + 1) . (48) 

Interestingly there exists a narrow window of parametric value for p, in which on for % = 1,2,3 
are positive and 04 and 05 are negative. An illustrative plot is shown in figure 4 for a particular 
value of p. 



10 



0' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

10 20 30 40 50 

Figure 4: The functions, t ai (rj) f(r), for i = 1, 2, • • • ,5 are plotted as functions of r/ for p = 0.52. 
This corresponding values for a^s are : a\ = «2 = 03 = 0.52 , 04 = —0.15 and 05 = —0.41. The 
isotropic expansion corresponding to a±, 02, 03 is plotted in red. The contraction corresponding 
to «4, 0:5 are plotted in green and black respectively. 

To conclude, we have presented a class of brane configurations which shows BKL oscillations 
near their cosmological singularities. It will be worthwhile to look for the signatures of these 
oscillations in their dual descriptions. We hope to report on this in the near future. 
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